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Statistical  analysis  of  measured  polarimetric  clutter 
data  at  different  range  resolutions 


M.  Greco,  F.  Gini  and  M.  Rangaswamy 

Abstract:  This  paper  deals  with  the  statistical  modelling  of  radar  backscattering  from  sea  surface  at 
low-grazing  angles  in  high  resolution  radar  systems.  High-resolution  polarimetric  data  at  different 
range  resolutions  (60,  30,  1 5,  9  and  3  m)  are  analysed  to  highlight  the  differences  in  clutter  statisti¬ 
cal  behaviour  due  to  changes  of  resolution  and/or  polari.sation.  The  clutter  data  were  recorded  by 
the  I  FIX  radar  of  McMastcr  University  in  Grimsby,  Ontario,  Canada. 


1  Radar  and  data  description 

The  clutter  data  processed  in  this  paper  were  collected  using 
the  McMa.ster  University  IPIX  radar.  IPIX  is  an  X-band 
search  radar,  capable  of  dual-polarised  and  frequency-agile 
operation  [I].  The  characteristic  features  of  the  IPIX  radar 
are  summarised  in  Table  I. 

The  radar  site  is  located  to  the  cast  of ‘Place  Polonaise’  at 
Grim.sby,  Ontario  (latitude  43.2II4N,  longitude 
79.5985  W),  looking  at  lake  Ontario  from  a  height  of 
20  m.  The  data  of  the  OHGR  database  are  stored  in  222 
files,  as  I  Obit  integers.  There  are  like-polarisations.  HH 
and  VV  (Lpol),  and  cross-polarisations,  HV  and  VH 
(Xpol),  with  coherent  processing  being  received,  leading 
to  a  quadruplet  of  I  and  Q  values  for  Lpol  and  Xpol. 

Wc  analysed  many  files  with  different  range  re.solutions 
and  recorded  on  different  days.  In  this  work,  we  summarise 
the  results  relating  to  only  five  files,  containing  more  than 
8  million  samples,  recorded  on  4  February  1 998  at  about 
22:30  hours  (local  time).  They  are  reprc.scntative  of  most 
of  the  results  obtained  from  all  the  processed  data. 
Unfortunately,  there  is  no  available  information  about 
the  wind  and  wave  observations  for  these  data  sets.  The 
relevant  parameters  are  summarised  in  Fig.  I  and  Table  2. 

It  is  important  to  observe  that  in  each  file  the  range  rc.sol- 
ution  is  different  and  that  at  30  m,  the  illuminated  zone  is 
different  with  respect  to  that  illuminated  at  60.  1 5,  9  and 
3  in  range  rc.solution.  These  files  were  chosen  to  highlight 
the  differences  due  to  the  change  in  the  resolution. 

The  IPIX  receiver  has  two  operational  modes  depending 
upon  the  .selected  RF  pul.se  width  (PW).  When  the 
•system  is  operating  with  PW  >  200  ns,  a  5  MHz  filter  is 
used  to  limit  the  receiver  bandwidth  to  ~5  MHz.  When  at 
PW  <  200  ns,  this  filter  is  bypassed,  so  the  bandwidth  of 
the  receiver  is  ~50MHz  to  match  the  minimum  20  ns 
PW.  Therefore  for  data  collected  with  PW  <  200  ns,  the 
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noise  floor  is  about  lOdB  higher  than  for  data  collected 
with  PW  >  200  ns  [2]. 

2  Statistical  analysis 

2. 1  Statistical  models  of  clutter  amplitude 

Many  distributions  have  been  proposed  in  the  literature  to 
model  the  amplitude  probability  density  function  (PDF)  of 
high  resolutipn  non-Gaussian  clutter  ([3  14],  and  references 
therein).  Here,  we  compare  the  empirical  PDF  with  log-normal 
(LN).  Wcibull  (W),  K  and  generalised  K  (GK)  PDFs.  The 
analytical  expressions  for  these  PDFs  and  their  moments 
;»;«(/?)  =  E[r"]  are  reported  below  [7],  where  r  =  |r|  denotes 
the  clutter  amplitude  and  ;  its  complex  envelope. 


LN  model 


Moments 


.  Cl  n' 

m^{n)  =  S  exp— — 


(1) 

(2) 


where  <t  is  the  .shape  parameter,  3  the  scale  parameter  and 
u{r)  the  unit  step  function. 

lyeihull  model  (WJ 


PDF  Pr(/-)  =  ^/‘''cxp[-Q  ju(/)  (3) 

Moments  m/fin)  =  h'T^-  +  I )  (4) 

where  c  is  the  shape  parameter  and  h  the  scale  parameter. 
The  Rayleigh  PDF  is  a  particular  case  of  the  Wcibull  PDF 
for  £■  =  2. 
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Table  1;  Characteristics  of  the  IPIX  radar 


Transmitter 

Receiver 

Parabolic  dish  antenna 

TWT  peak  power:  8  KW 

Two  receivers 

Diameter:  2.4  m 

Dual  frequency  simultaneous 

transmission:  8.9-9.4  GHz 

Outputs:  linear,  1  and  Q 

Pencil  beam  width  (azimuth  resolution):  1.1 

H-V  polarisation,  agile 

Receiving  polarisations:  H-V 

Antenna  gain:  45.7  dB 

PW:  20-5000  ns  (real)  5000  ns 

(expanded)  32  ns  (compressed) 

Data  acquisition:  sample  rate  from  0  to 

50  MHz  outputs:  linear,  1  and  Q 

Cross-polarisation  isolation:  30  dB 

PRF:  from  0  to  20  KHz 

Quantisation:  10  bit  up  to  16  bit 
effective  with  H/W  decimation 

Double  polarisation  with  central  feeder 

mtWifM  22J75J  AMST^J*  i  I9WII2M  22(NUV  ANTSTKt*  i  IMM2U4  22.>22l  AMVTF^ 


Numt  of  th« 
data  set 

Date,  time  of  ’  (iO/tM/IWK  22:37:53 
ac<j[uisition 

No.  of  ratine  cells  2K 

Mart  nyije  _ _ _ 

Kance  Ml  m 

resolution 

P\V _  _  _ _ 

Total  no.  of 

**!*?P.  _ _ _ _ 

Sample  for  MlllOO 

cell  santpled  at  Ml  m 

m 

Frv(|ueney 
RF 

Radar  and 
wave 
)>eomelry 


9.3V  OH/ 
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Fig.  1  Characteristics  of  the  anatyseci  files 
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where  !'(•)  is  the  gamma  function,  K,,  if  )  the  modified 
Bessel  function  of  the  third  kind  of  order  v  —  I,  u  the 
shape  parameter  and  /x  the  scale  parameter. 


GK  mode!  with  log-normal  texture 

r  r  2 

1  r  /  T  \i 

PDF  /Mr)=  - - 

•dllTrr  Jo  7^  L  7^ 

dr 

(7) 

Moments  m^(n)  =  (2m)"^~r^l  -f 

/  2 

(8) 

GK  model  with  f’eneralised  gamma  texture 

PDF  = 

..  .  .X  /■/i'\''/2^(i/  +  «/2f>)l'(l +n/2) 

Moments  w„(n)=(  — )  V - = - 

Vv/  r(u) 


dr 

(9) 

(10) 


Apart  from  LN  model,  the  other  PDFs  considered  here 
belong  to  the  so-called  compound-Gaussian  family. 
Specifically,  they  arise  from  the  product  of  two  random 
variables:  the  texture  and  the  speckle.  According  to  this 
model,  the  samples  of  the  complex  envelope  of  the  sea 
clutter  process  can  be  represented  as 


where  a  is  the  shape  parameter  and  m  the  scale  parameter. 


-(»■)  =  y/TU)x[i).  /■  =  1 , 2 . A/s  (II) 


Table  2:  Estimated  parameters,  60  m 


Cell 

W 

LN 

K 

LNT 

GK 

C 

b 

f> 

h 

V 

M 

m 

V 

A 

b 

VV-15th 

1.391 

0.014 

0.652 

0.010 

1.331 

1.268  X 

10 

4 

0.736 

8.760  X 

10 

5 

5.500 

2.798  X 

10 

$ 

0.461 

Mean 

1.293 

0.021 

0.695 

0.015 

1.119 

4.718  X 

10 

4 

1.020 

2.694  X 

10 

4 

20.57 

2.604  X 

10 

4 

0.427 

HH-15th 

1.324 

0.016 

0.677 

0.012 

1.122 

1.671  X 

10 

4 

0.867 

1.083  X 

10 

4 

6.010 

1.888  X 

10 

& 

0.408 

Mean 

1.226 

0.022 

0.722 

0.016 

0.927 

4.162  X 

10 

4 

1.133 

2.210  X 

10 

4 

15.55 

2.758  X 

10 

4 

0.428 

HV-15th 

1.386 

0.014 

0.654 

0.011 

1.316 

1.321  X 

10 

4 

0.744 

9.106  X 

10 

4 

5.341 

3.123  X 

10 

S 

0.465 

Mean 

1.292 

0.021 

0.695 

0.016 

1.117 

4.179  X 

10 

4 

0.980 

2.330  X 

10 

4 

19.70 

2.734  X 

10 

4 

0.430 

LNT,  LN  texture 
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w  here  i\  denotes  the  number  of  samples.  The  term  .v(/)  = 
V|(/) +7.Yy(/)  represents  a  stationary  complex-Gaussian 
process,  usually  called  speckle,  which  accounts  for  local 
back.scattering;  X]{i)  and  xqU)  are  the  in-phase  and  quadra¬ 
ture  components  of  the  speckle  complex  envelope  .y(/). 
Thc^'  satisfy  the  property  E{x,{i)]  =  fUplOl  =  0,  and 
^l'T(/)l  =  fl-Vpl/))  =  1/2,  so  that  £'(W;)|‘|  =  I,  that  is, 
the  speckle  complex  samples  have  unit  power.  The  factor 
Tii)  is  a  non-negative  real  random  process,  usually  called 
texture,  that  models  the  local  clutter  power. 

Given  the  speckle  and  texture  PDFs,  it  is  possible  to 
obtain  the  amplitude  PDF  as  follows 


r(/)  =  v/7^v(/)=>r(/)  =  |z(/)|  =  v/^W/)| 


-i; 


T)/j,(r)dT 


(12) 

(13) 


Fig.  2  Chiller  ampliliide  PDF,  W  polarisalion,  I5lh  range  cell, 

61)  III 


Fig.  3  NormaU.ied  duller  momenis,  FF  polarisalion,  15th  range 
cell,  61)  III 


The  parameters  of  the  theoretical  PDFs  are  estimated  by 
the  classical  method  of  moments  (MoM)  [15],  which  con- 
si.sts  of  equating  the  first-  and  second-order  experimental 
moments  with  the  corresponding  theoretical  moments. 
The  estimated  moments  are  given  by 

I  , 

"V«)=^l]k(/)r  (14) 

"s  ,=  | 

We  processed  =  60  000  samples  for  each  range  cell.  As 
concerning  the  parameters  of  the  GK-PDF,  we  encountered 
a  numerical  problem  with  the  above  approach.  Therefore  in 
this  case,  we  used  a  number  of  empirical  moments  greater 
than  the  number  of  unknowns.  As  a  result,  the  parameters 
V  and  b  were  c.stimated  as 


(v.  h)  =  arg  min7(v,  h)  =  arg  min  ^ 

iv,h)  I,,h)  m^in) 

(15) 

The  ab.solute  minimum  of  the  functional  in  (15)  was 
found  by  two  successive  two-dimensional  grid  searches, 
as  de.scribed  in  Farina  el  al,  [7].  Once  v  and  h  have  been 
estimated,  fi  is  obtained  from  an  estimate.  m/,(  1 ),  of  the 
first-order  moment. 

2.2  Range  resolutions  of  60,  30  and  15  m 

The  results  of  the  stati.stical  analysis  by  means  of  histograms 
and  moments  calculation  reveal  that  the  GK-PDF  yields  a 
good-fit  for  both  like-  and  cross-polarised  data  and  for  all 
the  three  resolutions.  Therefore  the  analysed  clutter  process 
can  be  accurately  modelled  by  a  compound-Gaussian 
process  with  GK-PDF,  provided  that  the  size  of  the  range 
resolution  cell  is  >  15  m  (note  that  the  Gaussian  model  is  a 
particular  case  of  the  GK  model). 

In  Figs.  2  and  3,  we  report  the  histogram  and  the 
moments  for  the  15th  range  cell,  VV  data,  60  m  range  res¬ 
olution.  The  numerical  results  for  the  other  range  cells 
and  the  two  other  range  resolutions  are  very  similar. 
Therefore  they  are  not  reported  here. 

In  Table  2.  we  report  the  mean  values  of  the  parameters 
estimated  for  each  theoretical  PDF.  The  results  show  that, 
for  a  resolution  of  60  m,  on  average,  the  HU  component  is 
spikier  (c  =  1.226)  than  both  VV  (f  =  1.293)  and  VH 
(c  =  1 .292)  components.  (Results  for  different  polari.sations 
have  been  compared  with  respect  to  the  mean  value  (c)  of 
the  estimates  c  of  parameter  c  of  the  Wcibull  distribution,* 
because  the  meaning  of  this  parameter  is  quite  easy  to 
understand.)  Moreover,  the  parameters  estimated  when  the 
range  resolution  is  30m  (Table  3)  show  that  the  data  are 
spikier  at  30  m  range  resolution  than  at  60  m;  this  was 
found  for  all  polarisations.  Moreover,  we  found  that  VV 
data  (c  =  1.094)  and  VH  data  (c  =  1.093)  are  .spikier  than 
HH  data  (f  =  1.218). 


Table  3:  Estimated  parameters,  30  m 


Mean 

w 

c 

b 

LN 

(T 

6 

K 

V 

A 

LNT 

.y* 

m 

GK 

V 

A 

b 

W 

1.094 

0.011 

0.681 

2.331  X  10  * 

1.573 

1.033  X  10  " 

12.80 

1.946  X  10  " 

0.436 

HH 

1.218 

0.003 

0.944 

1.701  X  10  = 

1.182 

8.966  X  10  ® 

29.77 

8.766  X  10  ® 

0.335 

HV 

1.093 

0.016 

0.011 

0.678 

2.468  X  10  * 

1.577 

1.094  X  10  * 

12.28 

2.061  X  10  ‘ 

0.437 

LNT,  LN  texture 
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Table  4:  Estimated  parameters,  15  m 


Mean 

w 

c 

h 

LN 

(T 

s 

K 

V 

A 

LNT 

m 

GK 

V 

A 

h 

W 

0.933 

0.029 

0.021 

2.309 

3.450  X  10  * 

3.602 

0.001 

IB 

HH 

0.874 

0.024 

2.592 

2.305  X  10  ^ 

3.281 

9.979  X  10  * 

HV 

1.018 

0.011 

1.905 

5.198  X  10  ® 

3.845 

1.078  X  10  * 

111 

LNT,  LN  texture 


Fig.  4  Chiller  ampliluiie  PDF.  FF  anil  HH  polarisalion\.  jiflh 
i  iinge  cell,  »i 


For  the  range  resolution  of  1 5  m  and  for  co-polarisations, 
we  observed  the  presence  of  some  more  heavy-tailed 
behaviour  (c.g.  in  the  7th  and  24th  cells  of  the  file 
1 9980204_223220),  but  generally  the  GK  model  provides 
a  good-fit  to  the  data.  The  values  of  the  estimated 
parameters  confirm  that  the  clutter  becomes  spikier  when 
range  resolution  increases  (i.e.  the  size  of  the  resolution 
cell  decreases).  We  also  noticed  that  on  an  average  HH 


data  are  spikier  (c  =  0.874)  than  VV  data  (c  =  0.933)  and 
VH  data  (c  =  1 .0 1 8);  the  same  happened  for  the  60  m  resol¬ 
ution  data  (Table  4). 

2.3  Range  resolutions  of  9  and  3  m 

Examining  the  histograms  obtained  by  analysing  the  file  at  a 
resolution  of  9  m,  wc  found  that  many  cells  of  co-polari.scd 
data  exhibit  heavy  tails  and  none  of  the  propo.scd  models 
yields  a  good-fit  to  the  data.  One  of  thc.se  cells  is  the  fifth, 
plotted  in  Fig.  4  for  VV  and  HH  data.  On  the  contrary, 
we  observed  that  for  cross-polarisations,  the  clutter 
process  can  still  be  accurately  modelled  by  a 
compound-Gaussian  process  with  GK-PDF.  Again,  the 
values  of  the  estimated  parameters  show  that  HH  data 
have  the  spikiest  behaviour  (HH:  c-  =  0.99l,  VV: 
c  =  1. 099,  VH:  c=  1. 1 75)  (Table  5).  The  same  was 
found  for  both  60  and  I5m  range  resolution  data. 

The  results  obtained  for  co-polarisations  at  a  range  resol¬ 
ution  of  3  m  do  not  show  significant  differences  with  respeet 
to  the  results  obtained  at  9  m.  Conversely,  the  analysis  for 
VH  polarisation  presents  some  dilTerence.  There  are  eells 
showing  histograms  with  tails  longer  than  the  average 
length  reeorded  at  lower  resolutions;  in  this  ease,  the 
eompound  model  eannot  be  used  to  model  elutter  data. 

The  estimates  of  the  parameters  are  reported  in  Table  6. 
The  results  show  again  that  HH  data  (c- =  1.307)  are 
spikier  than  VV  (c  =  1.542)  data.  With  respect  to  the 
other  resolutions,  estimated  values  of  c  arc  slightly  higher, 
probably  because  of  thermal  noise  elTeet, 


Table  5:  Estimated  parameters,  9  m 


Cell 

W 

LN 

K 

LNT 

GK 

c 

h 

fr 

S 

V 

A 

1-2 

m 

V 

A 

b 

VV-fifth 

0.585 

0.021 

1.207 

0.016 

0.129 

0.002 

4.861 

2.032  X  10  * 

0.216 

0.005 

0.638 

Mean 

1.099 

0.028 

0.806 

0.020 

0.781 

8.523  X  10  ‘ 

1.732 

3.261  X  10  ‘ 

35.69 

4.025  X  10  * 

0.2^7 

HH-fifth 

0.548 

0.015 

1.261 

0.011 

0.108 

0.002 

5.395 

1.041  X  10  ‘ 

0.147 

0.003 

Mean 

0.991 

0.021 

0.868 

0.015 

0.564 

5.852  X  10  * 

2.150 

1.825  X  10  ^ 

22.59 

3.777  X  10  ‘ 

HV-fifth 

0.626 

0.007 

1.153 

0.005 

0.154 

2.085  X  10  ‘ 

4.352 

2.367  X  10  ® 

0.459 

2.680  X  10  * 

Mean 

1.175 

0.009 

0.765 

0.007 

0.957 

8.157  X  10  ® 

1.461 

3.578  X  10  ® 

33.01 

3.832  X  10  ® 

0.259 

LNT,  LN  texture 

Table  6: 

Estimated  parameters. 

3  m 

Mean 

w 

LN 

K 

LNT 

GK 

C 

b 

(> 

s 

V 

A 

,y* 

m 

V 

A 

b 

W 

BSI 

B 

0.668 

0.028 

2.621 

0.001 

0.889 

6.470  X  10  * 

59.17 

1.997  X  10  ® 

0.148 

HH 

mm 

mm 

0.709 

0.019 

1.588 

5.580  X  10  ‘ 

1.118 

3.071  X  10  * 

54.83 

2.672  X  10  “ 

0.139 

HV. 

WM 

BB 

3.960 

4.550  X  10  ® 

0.608 

3.284  X  10  ® 

63.98 

2.01  X  10 

0.148 

LNT,  LN  texture 
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Table  7:  Probability  of  Type  I  error,  KS  test,  seventh 
range  cell,  VV  data,  9  m 


w 

LN 

K 

GK 

LNT 

52% 

0% 

51% 

33% 

0% 

LNT,  LN  texture 


The  .stati.stical  analysis  of  clutter  amplitude  was  com¬ 
pleted  by  applying  the  Kolmogorov-Smirnoff  (KS) 
goodness-ot-fit  test.  This  test  is  largely  used  to  determine 
which  distribution  provides  the  be.st-fit  to  the  data. 
Unfortunately,  in  some  cases,  it  was  not  useful  in  dis¬ 
tinguishing  between  different  long-tailed  models, 
bccau.se  it  places  an  equal  importance  on  all  regions  in 
the  probability  space.  Therefore  in  the  heaviest  part  of 
the  PDF,  that  is,  the  most  affecting  the  results  of  the 
test  is  the  ■bell'  area  or  the  body  of  the  PDFs,  many 
of  the  tested  PDFs  are  very  similar.  The  test  is  charactcr- 
i.sed  by  the  probability  of  Type  I  error.  The  Type  I  error 
is  the  probability  of  observing  under  H,)  a  sample 
outcome  at  least  as  extreme  as  the  one  observed  and 
hence  provides  the  smallest  level  at  which  the  observed 
sample  statistic  is  significant  [16]. 

For  resolutions  of  60  and  30  m,  for  all  the  tested 
models,  we  obtained  a  probability  of  Type  I  in  the 
range  of  I  5%.  Some  differences  were  found  at  higher 
resolutions.  In  Fig.  5,  we  show  the  empirical  and  theor¬ 
etical  distributions  for  a  heavy-tailed  range  cell,  the 
seventh  cell  for  VV  data  at  a  resolution  of  15  m.  The 
results  of  the  KS  test  for  this  cell  are  summari.scd  in 
Table  7.  For  resolution  higher  than  30  m,  the  KS  test  pro¬ 
vides  the  highest  probability  of  Type  I  error  particularly 
for  the  LN  distribution. 


3  Cumulant  domain  analysis 

To  perform  additional  analysis  of  the  compound-Gaussian 
model  and  to  investigate  whether  the  deviation  from  the 
theoretical  models  in  the  highest  two  resolutions,  that  is,  9 
and  3  m,  may  be  due  to  the  presence  of  non-ncgligible 
thermal  noise,  we  applied  the  theory  of  cumulants  [17].  It 
is  widely  known  in  the  literature  that  cumulants  of  order 
greater  than  two  for  a  Gaussian  process  are  identically 
zero  [17,  18].  Thus,  if  we  consider  the  clutter  process 
:(i)  =  }ii)  +  v(i).  where  v{i)  is  a  Gaussian  process  and 


i’(;)  is  a  non-Gaussian  process,  independent  of  v(i).  we  have 
. ^  =  . +  . 4  i) 

,  ( 1  o ) 

=  q(/| . 4_|).  for^  >2 

so  the  cumulants  of  if/)  can  be  derived  from  the  cumulants 
of  :{i).  In  our  case,  the  in-phasc  (1)  and  quadrature  (Q) 
components  of  the  thermal  noise  arc  zero-mcan-Gaussian 
processes,  then  only  non-Gaussian  clutter  contributes  to 
the  third-,  fourth-  and  fifth-order  cumulants  of  the  observed 
complex  data. 

We  estimated  from  the  data  the  second,  third  and  fifth- 
order  cumulants  at  zero  lags,  that  is,  for  /|  =  /:  =  /<  i  =  0. 
Then,  the  cumulants  are  normalised  with  respect  to  the 
second  order  cumulant  as  follows 

cl(0.0 . 0)  q^(0.() . 0) 

*  (4(0))*/-  (£V(0))*/= 

where  superscripts  1  and  Q  refer  to  the  in-phasc  and  quadra¬ 
ture  components,  that  is,  the  real  and  imaginary  parts  of  the 
complex  data.  We  compared  the  estimates  with  the  (normal¬ 
ised)  theoretical  cumulants  of  the  compound-Gaussian  model 
calculated  at  zero  lags.  All  the  theoretical  cumulants  of  odd 
order  calculated  at  the  origin  arc  equal  to  zero. 

In  Figs.  6  and  7,  we  show  the  normalised  cumulants  fiy 
and  fif  against  the  second-order  cumulant  for  the 


Hsiiinuleii  2nd  order  cumuluni 


Hsiimuled  2nd  order  cumulunt 


Fig.  6  Normalised  third-  ami  fifth-order  cumulants  against 
second-order  cumulant,  9  m 
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resolutions  of  9  and  3  m,  respectively.  The  results  show  that 
at  ii  range  resolution  of  9tti,  the  compound-Gaussian  model 
is  still  accurate.  Moreover,  it  can  be  adopted  to  mrxiel  the 
clutter  process  in  most  cells  also  for  VH  polarisation.  In 
fact,  for  VH  data  /aj  and  /r,  arc  close  to  zero.  Conversely, 
the  fifth-order  cumulant  shows  a  large  deviation  from’ 
zero  in  most  cells,  for  both  HH  and  VV  polarisations. 
At  a  range  resolution  of  3  m,  for  most  cells  and  all  polaris¬ 
ations,  the  estimated  cumulants  deviate  significantly  from 
zero.  This  is  an  indication  that  the  thermal  noise  is  not  the 
cause  of  the  deviation  from  the  compound-Gaussian  family. 


4  Correlation  analysis  and  power  spectrum 
estimation 


The  compound-Gaussian  clutter  model  assumes  the  pre¬ 
sence  of  two  components,  speckle  and  texture,  with  very 
different  correlation  times  (some  milliseconds  for  the  first 
component  and  some  seconds  for  the  .second  one).  If  the 
two  components  are  statistically  independent,  the  overall 
autocorrelation  function  is  the  product  of  the  autocorrela¬ 
tion  functions  of  the  two  components  [7,  19] 

/?/(/»)  =  £{r(/)r’(/ -I-  m))  = 

=  2Ry-^{m)  (/? V, (w)  -l-y/?,,, ,Vy('«))  ( 1 8) 

where  wc  exploited  the  fact  that  if  the  speckle  is  a  complex¬ 
valued  stationary  circular  process,  then  Rx,(ni)  =  R^  (m).  In 
practice,  the  decorrelation  time  of  the  signal  :U)’  called 
coherent  signal,  is  equal  to  that  of  the  fa,ster  component  [7], 


4. 7  Estimation  of  the  speckle  autocorrelation 
and  cross-correlation  sequences 

As  the  texture  can  be  considered  constant  over  short  time 
intervals,  we  can  estimate  the  speckle  autocorrelation  func¬ 
tions  7? v,(nt)  and  /?.v,a,/w)  by  using  coherent  signal  samples 
from  such  short  intervals  with  or  without  overlapping 


1 

‘A  1 

II 

L*=i 

<=« 

1] 

-  I 

1 

v.  l-m 

^  =  TT 

N^, 

E  -*(')-*('+"') 

/=(i 

] 

where  Ah  is  the  number  of  data  bursts  and  f*  the  estimated 
value  of  the  texture  in  the  itth  burst 


(21) 

\vhcrc  c*(0  =  -((A—  I  )A  -f  /).  We  used  three  different  values 
of  N:  64,  128  and  256,  with  an  overlap  between  bursts  of 
50%.  We  did  not  find  apparent  differences  in  the  estimation 
for  ^  =  64,  128  and  256,  so  we  show  here  only  the  results 
for  A  =128. 

Figs.  8  and  9  show  two  plots  of  the  real  and  the  imaginary 
parts  of  the  speckle  autocorrelation  function.  The  whole 
results  provide  a  clear  indication  that  for  all  resolutions 
and  all  polarisations,  the  speckle  correlation  time  is  about 
10  ms  long  and  the  behaviour  is  oscillatory,  as  in  Figs.  8 
and  9  for  60  and  3  m.  re.spectivciy. 


Fig.  8  Speckle  correlalioii  coefficients.  HU  polarisation,  first 
cell.  60  in 


4.2  Estimation  of  the  texture  autocorrelation 
sequence 

To  check  the  validity  of  the  hypothesis  made  over  the  cor¬ 
relation  times  of  the  two  components,  we  estimated  the 
texture  autocorrelation  sequence  with  the  formula 


It  is  useful  to  observe  that  with  a  50%  of  overlap,  we  can 
c.stimate  the  texture  correlation  and  covariance  every  A/2 
lags.  To  estimate  the  texture  correlation,  wc  .set  A=  128. 
In  the  figures,  we  plot  the  texture  correlation  coefficient 
that  is 


_  C^(A/w/2) 
2  /  C.(0) 


(24) 


Figs.  10  12  show  the  texture  correlation  coefficient.  At 
the  same  resolutions  and  without  differences  in  the  polarisa¬ 
tions,  the  texture  correlation  time  is  on  the  order  of  .seconds. 
Furthermore,  the  texture  presents  periodicities  with  a  period 
of  8  s  at  a  range  resolution  of  60  m  and  of  3  s  at  a  range  ros- 
olution  of  30  m.  The  periodicity  is  particularly  evident  in  the 


Fig.  9  Speckle  correlation  coefficients.  VV  polarisation  eivhth 
cell.  3  m 
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Fig.  10  Texture  miUK  Ovariance  funclion.  first  cell  N=  I2K.  60  m  Fig.  12  Texture  autueuvariuiiee  function.  i'F  poluri.sation. 

eighth  cell,  .1  in 


VV  polarised  data  for  the  resolution  of  60  m  in  all  the  files 
we  analysed. 

Our  results  show  that,  with  increasing  resolution,  the 
texture  correlation  time  decreases,  but  still  in  the  order  of 
a  few  seconds  and  the  periodicities  tend  to  disappear,  due 
to  the  -Strong  contribution  of  the  thermal  noise  (Figs.  9 
and  10). 

Figs.  1.1  and  14  report  an  example  of  the  average  spectro¬ 
gram  in  semi-logarithm  scale  for  each  considered  cell, 
calculated  as 


.r(/)e-'*"  . 

.  +  > 

K  (25) 

where  is  the  number  of  sequences  in  which  the  received 
vector  for  each  cell  has  been  divided,  the  number  of 
samples  per  sequence,  k  the  normalised  frequency  and 
Piik)  is  the  Ath  sample  of  the  periodogram  of  the  /1h 
sequence. 

The  periodogram  shows,  for  all  the  polarisations,  for  all 
the  range  resolutions,  a  peak  located  around  I  SO  Hz. 
Moreover,  with  a  resolution  of  60  and  30  m,  all  the  analysed 
cells  show  a  bimodal  spectrum,  particularly  evident  in  HH 
polarisation,  and  then  a  second  peak  near  —  1 50  Hz;  the 
power  of  the  .second  peak  is  much  lower  than  the  power 
of  the  main  one  (Fig.  I3).  From  the  resolution  of  ISm, 
the  IPIX  radar  seems  to  add  a  frequency  interfering  line 
in  the  spectrum  at  about  —220 Hz  (Fig.  I4).  The  line  at 


,V  .V 

I  1  1 

r=l  ^  r=l  '  c 


rN, 

E 

Mr-IW, 


A  =  1.2. 


0  Hz  is  due  to  a  residual  of  the  continuous  component.  It 
is  evident  that,  as  resolution  increases,  the  thermal  noise 
effect  becomes  very  important,  in  fact,  from  our  results, 
only  partially  shown  here,  it  is  evident  that  the 
cluttcr-to-noise  ratio  (CNR)  decreases  from  ~24dB  for  a 
resolution  of  60  m  with  VV  data  to  less  than  -  5  dB  for  a 
resolution  of  3  m  with  VH  data.  The  CNR  has  been 
roughly  estimated  from  the  spectrum  figures  reading  the 
value  of  the  noise  floor  from  each  figure  and  calculating 
clutter  power  as  the  difference  between  the  overall 


Fig.  13  Average  /tower  spectra!  density,  first  cell.  HH  /tolaris- 
alion.  60  in.  CNR  ^  24  dB 


Fig.  11  Texture  autocovariance  function.  HH  /tolarisation.  first 
ccli.  Ifi  III 


Fig.  14  Average  power  spectra!  densitv.  VV  polarisation,  eighth 
cell.  .1  m.  CNR  ^  4  dB 
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Fig.  15  Mi'an  rungi-  lexliirc  aiilocovariance.  ami  HH 
polarisulion 

disturbance  power  and  noise  power.  The  values  calculated 
for  each  cell  is  reported  in  the  figure  captions. 

4.3  Mean  range  texture  autocovariance  sequence 

To  conclude  the  correlation  analysis,  in  order  to  highlight 
further  differences  due  to  the  resolution,  we  also  calculated 
the  average  range  autocovariance  function  of  the  texture 
given  by 

I  . 

'  ^  m=  I 

1  ^'b 

=  VaT Z)  Z +  ”>  -  <26) 

‘  t>  l'  /tt=  I  ;■=:  I 

where  is  the  estimate  of  the  texture  on  the  wth  burst  of 
the  /th  cell,  the  number  of  illuminated  cells,  Nt,  the 
number  of  bursts  and  t„,  the  texture  average  value  in  the 
/nth  burst 

I 

=  xrZ^'"<'>  <27) 

'’'c  (=1 

Because  at  30  m  range  resolution,  the  illuminated  zone  is 
different,  we  only  compared  the  results  found  at  15,  9  and 
3  m  range  resolutions. 

Figs.  10  and  1 5  show  the  results  obtained  for  VV  and  FIH 
polarisations.  From  the  figures  it  is  evident  that  with  a  res¬ 
olution  of  3  m,  it  is  possible  to  highlight  and  resolve  shorter 
range  periodicities  that  are  not  visible  in  the  other  resol¬ 
utions  in  all  the  polarisations. 

5  Conclusions 

Our  results,  based  on  the  analysis  of  many  files  at  different 
range  resolutions  of  the  Grimsby  data  set,  recorded  by  I  FIX 
radar,  reveal  that  the  GK-PDF  provides  a  good-fit  to  the  data 
for  all  polarisations  with  resolutions  of  60  and  30  m. 
Increasing  the  range  resolution  up  to  15  m.  the  model  con¬ 
sistently  yields  a  good-fit  in  most  cells  for  co-polarisations 
even  if  the  presence  of  longer-tailed  cells  with  respect  to  the 
lower  resolutions  has  been  observed.  The  model  is  valid  for 
all  VFI-FIV  polarised  data.  At  9  and  3  m,  both  VV  and  FlFl 
polarisation  histograms  exhibit  very  long  tails,  so  that  the 
compound-Gaussian  model  is  not  valid  any  longer. 

4S() 


Moreover,  at  these  resolutions,  the  thermal  noise  contri¬ 
bution  becomes  non-negligibic.  The  GK  model  provides  a 
good-fit  in  almo.st  all  cells  for  VH  polarisation  at  9  m,  but 
at  3  m  also  the  cros.s-polarisations  prescjii  long-tailed  cells 
and  the  compound  model  fails  to  provide  a  good-fit. 

Estimated  parameters  analysis  shows  that  HH  component 

IS  the  most  spiky  one  at  60,  I  5,  9  and  3  m  range  resolution; 

at  a  resolution  of  30  m,  on  average,  an  HH  component  is  less 
spiky  than  both  VV  and  VH  components  (remembering  that 
at  30  m  range  resolution  the  illuminated  zone  is  different). 

Cumulants-domain  analysis  confirms  the  results  obtained 
by  the  histograms  analysis.  At  9  m  range  resolution 
compound-Gaussian  model  provides  a  good-fit  in  most 
cells  for  VH  polarisation.  The  fifth-order  cumulant  shows 
a  large  deviation  from  zero  in  most  cells  for  both  HH  and 
VV  polarisations. 

At  3  m  range  resolution,  most  cells  in  all  polarisations 
present  estimated  values  of  the  cumulants  very  different 
from  zero.  In  this  instance,  the  deviation  from  the 
compound-Gaussian  model  is  not  due  to  the  presence  of 
thermal  noise. 

About  the  correlation  and  spectral  analysis,  based  on  our 
results,  we  can  conclude  that  in  all  the  analysed  resolutions, 
in  each  polarisation,  average  speckle  correlation  time  is 
~  10  ms.  Moreover,  without  differences  in  the  polarisations, 
the  texture  correlation  time  is  .some  seconds  long;  the 
texture  presents  periodicities  with  a  period  of  8  s  at  a 
range  resolution  of  60  m  and  of  3  s  at  a  range  resolution 
of  30  m.  The  periodicity  is  particularly  evident  in  the  VV 
polarised  data  for  the  resolution  of  60  m.  With  increasing 
resolution,  jhc  texture  correlation  time  becomes  shorter, 
but  still  on  the  order  of  a  few  seconds,  and  the  periodicities 
tend  to  disappear,  due  to  the  strong  contribution  of  the 
thermal  noise.  Differences  in  HH.  VV  and  VH  time  covari¬ 
ances  arc  not  evident.  Only  in  range,  the  VH  texture  exhibits 
a  longer  correlation  time  than  the  like-polarised  data. 

The  analysis  of  the  periodogram  shows  that  for  all  the 
polarisations  and  all  range  resolutions,  a  peak  is  located 
around  1 50  Hz.  Moreover,  as  the  resolution  increases,  the 
thermal  noise  effect  becomes  very  important;  the  CMR 
decreases  from  ^24  dB  for  a  resolution  of  60  m  and  VV 
data  to  less  than  —  5dB  for  a  resolution  of  3  m  and  VH 
data.  Then,  the  presence  of  the  thennal  noise  at  very  high 
resolution  strongly  affects  the  behaviour  of  the  recorded 
disturbance. 
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